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We investigate Choptuik scaling in the spherically symmetric collapse of a mass- 
less scalar field in higher dimensions using Painleve-Gullstrand (PG) coordinates. 
Our analysis gives reliable results for the critical exponent and period echoing up 
to seven dimensions and confirms the presence in higher dimensions of the cusps 
in the periodic scaling relationship recently observed in four dimensional collapse. 
In addition, our value for the critical exponent in seven dimensions is consistent 
with that obtained by Bland et al, who argued that the critical exponent increases 
monotonically with dimension to an asymptotic value of 1/2. 
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1 Introduction 

The formation of a singularity at the center of black holes signals the breakdown of 
classical general relativity. It is speculated that a quantum theory of gravity would 
resolve the singularity problem. Some popular theories of quantum gravity, such as 
string theory, depend upon the existence of more than four dimensions. It has even 
been suggested that microscopic higher dimensional black holes could be produced in 
the LHC, and experimental bounds on their production have appeared in the literature 
p] . It is therefore important to understand the details of black hole formation in higher 
dimensions. 

In the early 90 's Matthew Choptuik [2] numerically simulated the collapse of a spher- 
ically symmetric massless scalar field and found a simple universal scaling relationship 
between physical observables describing the black hole on formation and the parameters 
describing the initial scalar field distribution. Since then, this universal scaling has been 
confirmed in a large number of different settings. (For a review of Choptuik scaling see 
[3].) To be more precise, as a single parameter, A, for any family of initial data, is varied 
the initial data is separated into two classes. For A less than some critical value A*, the 
matter bounces from the origin and ultimately dissipates. For A > A*, an event horizon 
forms and the end state is a black hole. As A — > A* from above, the final spacetime 
approaches arbitrarily closely the critical black hole solution which in effect describes a 
zero mass black hole that exhibits a naked singularity at the center. Near criticality, any 
geometrical quantity, X, characterizing the evolution of the black hole at the time of 
formation, obeys a scaling law of the form: 

ln{X) = jln\A - A, | + f(ln\A - A*\), (1.1) 

where 7 is the critical exponent and / is a periodic function with period, T. 7 and T 
are universal in the sense that they are independent of the quantity being measured and 
also of the type of initial data. They do depend on the type of matter that is collapsing, 
however. The specific form of the oscillatiing function, /, is known not to be universal. 
The most common quantities used to test this relationship are the mass of the final black 
hole (or equivalently the radius of the event horizon) and the maximum value of the 
curvature in subcritical collapse. Previous calculations for the collapse of a spherically 
symmetric massless scalar field in four spacetime dimensions in Schwarzschild, double 
null and Painleve Gullstrand (PG) coordinates have all given 7 0.37 and T rj 4.6 

SI El El El El ED] . 

In Schwarzschild and double null coordinates the form of the periodic function, /, is 
well approximated by a small amplitude sine function. By contrast, in PG coordinates 
the periodic function in four dimensions showed distinctive cusps [6J. This difference in 
the form of / in PG coordinates is perhaps surprising, but not inconsistent. The PG 
calculation did not measure the final black hole mass after all matter had fallen through 
the horizon, which is independent of slicing and would exhibit the same behaviour in 
all coordinates. Instead, what was plotted was the radius of the apparent horizon on 
formation, a quantity that does depend on the slicing. It is nonetheless a geometrical 
quantity that exhibits Choptuik scaling, as shown in [Sj. What is less clear is whether the 
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cusp-like features of the scaling function is a peculiarity of PG coordinates in 4 dimensions 
or generic in some sense. 

Most previous calculations of Choptuik scaling have involved four dimensional black 
holes. The first higher dimensional analysis was done by Garfinkle and Duncan 
More recently, a program was initiated [12] whose purpose was to calculate the critical 
exponent and echoing period for spherical collapse in arbitrary spacetime dimensions 
using the so-called dilaton gravity formalism. Preliminary results were obtained in 5 
and 6 spacetime dimensions. Subsequently, accurate results were obtained up to D=14 
[7J [8] . These results provide strong evidence that the critical exponent was a monotonic 
function of D that converged asymptotically to 1/2. Also around that time, Oren and 
Sorkin[13] produced results for the critical exponent up to D=ll which they used to 
speculate that the critical exponent was not a monotonic function, but instead achieved 
a maximum near D=10. Moreover, it was argued on the basis of the behaviour of black 
string solutions that a critical dimension near D=10 might not be unexpected. It should 
be noted that both sets of calculations were done using double null coordinates, using 
different parametrizations of the fields. 

The higher dimensional calculations referred to above were done in double null co- 
ordinates, whereas the purpose of the present paper is to look at spherically symmetric 
scalar field collapse in higher dimensions using Painleve-Gullstrand (P-G) coordinates. It 
is perhaps worth highlighting the relative merits of the two different sets of coordinates. 
The main advantage of double null coordinates for studying critical behaviour has to do 
with the fine spatial resolution that is possible due to the convergence of outgoing null 
rays near horizon formation. However, this same convergence makes it difficult, if not 
impossible to allow numerical simulations to run all the way up to horizon formation (see 
numerics section of [9]). For this reason the mass/horizon radius of the black hole in null 
coordinates must be approximated by choosing in advance how close one wishes to get to 
horizon formation. This problem is also encountered in Schwarzschild coordinates. PG 
coordinates provide spatial slices that are regular across future horizons, so that the code 
can be run up to and even beyond horizon formation. Away from criticality, this allows 
one to map out the structure of the trapping horizon as done in [6]. Near criticality, we 
can precisely (to numerical accuracy) determine the radius of the horizon on formation. 
Note that this quantity is different from the mass of the final black hole. It marks the 
location that the P-G slicing intersects the trapping horizon and therefore verifies the 
scaling relationship in a different geometrical quantity therbye providing independent 
measures of the critical exponent and echoing scale. This also explains why a different 
oscillatory function is possible. The main disadvantage of P-G coordinates is that one 
does not get an automatic improvement in the spatial resolution from the convergence 
of null rays. Moreover, as discussed in [5] the time steps are limited by the null cone 
structure. 

As with previous calculations, the resulting numerics get more difficult as the number 
of spacetime dimensions is increased. We obtain reliable results in 4, 5, 6 and 7 dimensions 
for the scaling law obeyed by the areal radius of the apparent horizon on formation. We 
investigate the periodic function, /, in this scaling relation and confirm the cusp-like 
nature in higher dimensions. In addition, the higher dimensional critical exponents and 
periods are calculated and compared to previous work. Our result in 7 dimensions is 



2 Equations of Motion 



consistent within error to those of [7] and inconsistent with JT3] 

The paper is organized as follows. Section II reviews the formalism and derivation of 
the evolution equations in P-G coordinates for spherically symmetric collapse of a scalar 
field in arbitrary spacetime dimensions. Section III describes the numerical method and 
presents our results. Section IV closes with conclusions and prospects for future work. 



2 Equations of Motion 

The equations of motion are derived from the action for gravity minimally coupled to a 
massless scalar field, 



S\g D ,M = y^q^ J dx D ^g~D)R D (g D ) - ± J dx D v ^o)\ V^d| 2 , (2-1) 

where D is the number of dimensions and Gp, g^, Rd and ipo are the D dimensional 
gravitational constant, metric, Ricci scalar and the massless scalar field. |V ' d^d\ 2 is short 
hand for g^d^ D d v ^) D - 

2.1 Dilaton Formalism 

Since we are concerned with spherical symmetry, dimensional reduction and a field repa- 
rameterization can be used to put the action and equations in so called 2-D "dilaton" 
gravity form. This simplifies the problem and allows a unified treatment of all space- 
time dimensions. See [H] for a review and [15] for its applications to black holes. The 
dimensionally reduced action in dilaton form is: 

S[gAM = \ j dxdty/{^g)[(l/G){ ( f>R{g)+V{ ( f>))-h{ ( f>)\V^\ 2 ], (2-2) 
It is obtained by making the following field definitions: 



n n7T n x n xrfj.\ f H^l - 1) 

(n — \)v [n > 8(n — 1) \n(n — lj 

hty) = 8(n-1) 0, ^ = v / ^>)^ D , (2.3) 
n Jo 

where r is the areal radius, G and ip are the spherically symmetric versions of the grav- 
itational constant and the scalar field, n = D — 2 is the number of angular dimensions 
and the volume of a unit n-sphere. Note that the dilaton, 0, has a geometrical 

interpretation as the area of a sphere at fixed distance from the center. With these field 
redefinitions the physical metric is related to the two dimensional metric, ds 2 as follows: 

ds 2 

Using this formalism we can work with the 2D metric, g^ u {r,t), and associated Ricci 
scalar, R(r,t), without losing the effects of the D dimensional curvature. 
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2.2 Gauge Fixing 

The derivation of a partially gauge fixed hamiltonian from the action (I2.2p was carried 
out in [T7] using the general ADM metric, 

ds 2 := e 2p [-a 2 dt 2 + (dx + Ndt) 2 }, (2.5) 

where p, a and N are arbitrary functions of the spacetime coordinates, x and t. The 
gauge chosen is given by, 

j(<f>) = l<f>'. (2.6) 

This gauge choice forces the spatial coordinate to be the area radius, r, found in the 
Schwarzschild and PG metrics and gives the hamiltonian, 



dra 



+ dr 
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2p 
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(2.7) 



where II p and 11^ are the conjugate momenta corresponding to p and ip. The fields M. 
and Gm are given by, 



M 



1 

" 2G 
Gm = 



G 2 e~ 2 m 2 p - 



L/\2 



:,2p 



In the above M. is so-called the mass function that generalizes the Misner-Sharpe 
mass [TB] to higher dimensions. It aysmptotes to the ADM mass as r — > oo. Gm is the 
energy density of the free scalar field. 

We now completely fix the gauge using, 



(2.8) 
(2.9) 



(2.10) 



as was done in [6]. This gauge choice forces the metric into the non-static, higher dimen- 
sional generalization of the PG metric, given by 



ds 2 =m 



-a 2 dt 2 



dr 



I2GM 



crdt 



(2.11) 



An advantage of the PG coordinate system is that spatial slices t = constant are regular 
across apparent horizons that form during evolution. The Hamiltonian of Eq. fl2.7|) and 
gauge choices of equations (12.61) and (12.101) give the fully reduced equations of motion 



W = o \ 7 == — h 



(2.12) 
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-M 7 = Qm + 4>% 
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=<r = 0. 



y/2GMm 

Notice that these equations depend on the number of dimensions via Eq. (l2.3|) . 



(2.13) 

(2.14) 
(2.15) 



3 Methodology and Results 
3.1 Computational Details 

In order to solve the equations of motion we used the initial conditions: 



if; = Ar 2 exp 

n> = o, 



r — To 
B 



(3.1) 



where A, B and ro are parameters that can in principle be varied to study mass scaling. 
These initial conditions are the same as those used in [5J. We verified universality of our 
results by varying the amplitude, A, and the width, B, of the initial pulse. 

The equations of motion, (I2.12p -f !2.15j) . were solved using a forth order Runge-Kutta 
scheme with deriviatives being calculated using the finite difference method. 

In order to maintain stability we used an adaptive time step, At(t): 



At(t) = min r 



dt , 
— Ar r j 
dr 



(3.2) 



where Ar(r) is the spacing of the spatial lattice and ^ is the inverse of the local speed 
of an ingoing null geodesic. 

Much of the interesting behaviour of the collapse near criticality occurs near the origin. 
This requires close spacing of the spatial mesh near r = 0. It is not computationally 
realistic, however, to use this spacing along the entire spatial slice, which needs to be 
long enough so that none of the mass leaves during the simulation. The use of close 
spacing for the entire slice would dramatically increase the simulation time. For this 
reason the spacing near the origin was set to 10~ 5 and then smoothly increased to 10~ 2 
over the first 100 of the 1200 total lattice points. 
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3.2 Results 

To find Choptuik's mass scaling relation we first found the critical values, A* and B*, 
using a bisection method to an accuracy of A A* /A* m AB*/B* rs 10~ 15 . We then varied 
the parameters A and B in our initial data function, Eq.f l3.ip . and found the mass of 
the black hole at formation, Mbh, as a function of the initial parameters. An apparent 
horizon forms when an outgoing null geodesic becomes momentarily stationary. For a 
metric in P-G form as in Eq. fl2.lip . this is signalled by the condition: 

2GM-j^)\ rf = 0, (3.3) 

where rf is the areal radius of the horizon on formation. The mass of the black hole at 
horizon formation is then given by the value of the mass function -M(r/) at that point. 
We emphasize again that this is not the ADM mass of the final static black hole that is 
left behind once all the matter has fallen through the horizon. 

The mass scaling plots can be seen in Fig.(|T]). A straight line, which osculates the 
curve, has been plotted with the data to illustrate the linear term in the mass scaling 
relation, Eq. fll.il) . Notice that these plots confirm the existence in 5, 6 and 7 dimensions 
of the cusps in the periodic function that were originally noted in four dimensions in [5] . 
It is also important to note from Fig.fjT]) that as we move away from the near critical 
region (i.e. as A — A* gets large) the cusps systematically move away from the critical 
straight line. In 4 dimensions the results move above the critical line, in five dimensions 
they more or less stay on the straight line, but in higher dimensions the cusps move 
downward. This illustrates graphically that in higher dimensions the critical exponent 
will be underestimated if one is not sufficiently close to criticality. 

Points in the data that were close to the straight lines (< 2%) of Fig.fJT]) were chosen 
and were used to find the critical exponent using linear regression. An estimate of error 
was found by slightly varying which points were chosen to calculate the slope. The period 
of / in Eq.f ll.ip was calculated by measuring the distance between the cusps of the plots. 
A measure of error for these cusps was taken as the largest deviation from the average. 
All measurements were made using the first three periods of the data where the near 
criticality approximation is valid. This analysis was done for both of the cases where A 
and B were varied in the initial data. The results can be seen in table ([1]). These results 
show good agreement with the results of Bland et. al. [7] and Ziprick and Kunstatter 

In order to give a clearer sense of why our results are limited to seven dimensions, we 
present eight dimensional data in Fig.(j2]). In lower dimensions we were able to obtain 
supercritical data down to ln(A — A*) —17. In eight dimensions we were limited to 
ln(A — A*) < —12.6 without decreasing the lattice spacing near the origin, which in turn 
would dramatically increase the simulation time. The effect of this is seen in figure (J2J): we 
are unable to get close enough to criticality for the slope to be approximately constant 
over three periods. For this reason we were unable to reliably calculate the critical 
exponent for eight dimensions with the same accuracy as lower dimensional simulations. 
If we nonetheless obtain a slope from the first two periods of the simulation (as opposed 
to three) then we find 7 = 0.44 ± 0.02 and T = 3.0 ± 0.1. This result lies between that 
of Bland et. al. (7 = 0.4459 ± 0.0054, T — 3.11 ± 0.1) and Sorkin and Oren [T3] 
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(a) 4 Dimensions 





(b) 5 Dimensions 




(c) 6 Dimensions (d) 7 Dimensions 

Figure 1: Mass Scaling in 4, 5, 6 and 7 dimensions with the A parameter varied 



(7 = 0.436 ± 2%, T — 3.1 ± 3%). It should be noticed that given the arguments above 
and the negative curvature that we see in figure ([2]) our current value is almost certainly 
an underestimate. That is, a third cusp closer to criticality would most likely raise the 
value of the slope, bringing the critical exponent closer to that of [7j. 



4 Conclusion 

We have successfully confirmed the existence of cusps in the mass scaling function, 
Eq.f ll.il) . in 4 to 8 dimensions as observed in 4 dimensions by Ziprick and Kunstatter 
[6]. In addition the mass scaling plots obtained using PG coordinates give critical data 
which agrees with previous results in 4, 5 and 6 dimensions. However, in 7 dimensions, 
our critical exponent and echoing period agree with Bland et. al. [7] but disagree with 
the values obtained by Sorkin and Oren [15] . It is important to note that [7] claimed 



4 Conclusion 



Table 1: Critical Exponent, 7 and Period, T for Four to Seven Dimensions, D 



D 




7 


T 


4 


A varied 


0.378 ± 0.002 


4.4 ± 0.3 




B varied 


0.379 ± 0.002 


4.4 ± 0.2 




from [7] 


0.374 ± 0.002 


4.55 ± 0.1 




from [6] 


0.375 ± 0.004 


4.6 ± 0.1 




from [13] 


0.372 ± 1% 


4.53 ± 2% 


5 


A varied 


0.413 ± 0.002 


3.9 ± 0.8 




B varied 


0.416 ± 0.002 


3.7 ± 0.4 




from [7] 


0.412 ± 0.004 


3.76 ± 0.1 




from [13] 


0.408 ± 2% 


4.29 ± 2% 


6 


A varied 


0.429 ± 0.003 


3.4 ± 0.3 




B varied 


0.428 ± 0.002 


3.3 ± 0.2 




from [H] 


0.424 


3.03 




from 1/1 

11 Will 1 1 1 


430 4- 003 


3 47 + D 1 




from [13] 


0.422 ± 2% 


4.05 ± 2% 


7 


A varied 


0.440 ± 0.005 


3.1 ± 0.1 




B varied 


0.440 ± 0.006 


3.1 ± 0.4 




from [7] 


0.441 ± 0.004 


3.36 ± 0.1 




from [TBI 


0.429 ± 2% 


3.80 ± 2% 




-20 -15 


-10 


5 




-10- 



-12 - 



Figure 2: Mass Scaling in 8 dimensions (A varied) 

that the critical exponent was a monotonic function of spacetime dimension that asym- 
potes to 1/2, whereas [13] suggested that the critical exponent peaks near D=10. While 
both claims are intriguing, only one (at most) can be right. Since we are limited at the 
moment to D < 7 we cannot make any definitive claims about the asympotic behaviour. 
Our result in D = 7 supports that of [8], as does our more tentative result in D = 8. It 
is however impossible to make definitive claims about the asymptotic behaviour without 
pushing the PG calculation to higher dimensions. 
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